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We study the energy relaxation times (Ti) of superconducting transmon qubits in 3D cavities as a function of 
dielectric participation ratios of material surfaces. This surface participation ratio, representing the fraction of 
electric field energy stored in a dissipative surface layer, is computed by a two-step finite-element simulation 
and experimentally varied by qubit geometry. With a clean electromagnetic environment and suppressed 
non-equilibrium quasiparticle density, we find an approximately proportional relation between the transmon 
relaxation rates and surface participation ratios. These results suggest dielectric dissipation arising from 
material interfaces is the major limiting factor for the Ti of transmons in 3D cQED architecture. Our 
analysis also supports the notion of spatial discreteness of surface dielectric dissipation. 


Circuit quantum electrodynamics (cQED) systems 
have emerged as promising platforms for quantum infor¬ 
mation processing, powered by dramatic improvement of 
the coherence times of superconducting qubits over the 
past decade^. Such an improvement has been the result 
of collective efforts in multiple aspects^, such as suppres¬ 
sion of charge noise and flux noise^, better control of the 
electromagnetic environment^, elimination of deposited 
dielectric materials^, development in surface treatment^, 
dilution of surface effects by expanding field volume^, and 
improved filtering and shielding against stray radiation®. 
However, it has been difficult to quantify how much each 
of these individual measures contribute to the overall im¬ 
provement. As a result, it remains elusive what the domi¬ 
nant limiting factors are for the coherence of state-of-the- 
art superconducting qubits such as the 3D and planar 
transmons. 

The superior lifetimes (Ti) of qubits with larger 
footprints^ or with more advanced surface preparation® 
strongly suggest the important role of dielectric dissipa¬ 
tion^ from material surfaces. In this letter, we quantita¬ 
tively extract surface dielectric dissipation in transmon 
qubits through a combined experimental and numerical 
study. We find that surface dielectric dissipation is prob¬ 
ably still the major limiting factor for Ti of transmons in 
3D cQED architecture, and so far there is no indication 
of additional loss mechanisms (up to the level of Q ^ 10^) 
under our experimental condition. Our analysis also indi¬ 
cates that surface loss for a sub-micrometer area cannot 
be captured by a uniform loss tangent model, consistent 
with the hypothesis of discrete dissipation from a small 
number of microscopic two-level states (TLS)^“^^. 

Relaxation of superconducting qubits or resonators can 
be caused by many dissipative channels such as dielectric 
loss, conductive loss, and radiation into free space^. Di¬ 
electric loss can be further decomposed into contributions 
from various materials or components, so that: 
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where Ti , Q and oj are the relaxation time, quality factor 
(for energy decay) and angular frequency of the qubit 


or resonator, Tq is the relaxation rate induced by non¬ 
dielectric channels, Qi = 1/ian5i is the quality factor of 
the material with a dielectric constant of (with tan 5 
known as the loss tangent), and pi is its participation 
ratio defined as the fraction of electric field energy stored 
within the volume of this material. 

Crystalline substrates of cQED devices often store a 
large fraction of electric field energy {pi ~ 90%), but re¬ 
portedly show very small loss tangent (tanJ^ < 10“® for 
bulk sapphire^^ and silicon^). On the other hand, if a mi¬ 
croscopic layer of contaminants such as oxide, adsorbed 
water or organics forms at the metal-substrate (MA), 
substrate-air (SA) and metal-air (MA) interfaces^®’^^, 
they have much smaller pi but may still induce significant 
dissipation with a large tand^ on the order of 10“®-10“^. 
Previous studies^®^^® have found a positive correlation 
between the quality factors of planar resonators and their 
feature sizes which can be used to vary pi. However, a 
quantitative test of Eq. (1) has been challenging due to 
the presence of other energy relaxation channels (Tq) that 
have not been fully under control. 

Here we study the energy relaxation time, Ti, of trans¬ 
mon qubits as a function of surface dielectric participa¬ 
tion ratio. Pi- Strong suppression of radiation loss is 
achieved by implementing the 3D cQED architecture"^ 
where the 3D cavity enclosure provides a clean electro¬ 
magnetic environment free of spurious modes. The cav¬ 
ity Q and qubit-cavity detuning are sufficiently large to 
avoid any appreciable Purcell effect. Qubit relaxation 
due to non-equilibrium quasiparticles can be estimated 
and suppressed by monitoring and controlling quasipar¬ 
ticle decay time^®“^^. Furthermore, transmons are less 
sensitive to vortex ac loss than linear resonators because 
most inductive energy is stored in the Josephson junction 
rather than the electrodes subjected to vortex penetra¬ 
tion. Suppression of these relaxation channels allows us 
to vary the qubit geometry to change pi by more than an 
order of magnitude, making quantitative comparison of 
surface dielectric loss in different devices viable. 

Each qubit in this study is composed of a single 
Al/A10a;/Al Josephson junction and a pair of electrodes 
forming a shunting capacitor. We report Ti measured 
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FIG. 1. Geometry of four different designs of trans- 
mon qubits used in this study. Most dimensions of the 
electrodes for each design (A, B, C, D) are fixed and labeled 
in respective panels (a, b, c, d). An exception is that Design 
C has multiple variations with the dimensions g = w ranging 
from 1.5 fim to 30 fim. For convenience, we define “leads” 
as the portion of the electrodes with widths 1 ^m or smaller, 
which only appears in direct connection to the junction in 
all our designs, as shown in dark brown. The rest of the elec¬ 
trodes wider than 1 fj,m are called “pads” and shown in lighter 
grey, (e) Schematic of the standard 3D cQED setup. Trans- 
mon qubits are installed in rectangular waveguide cavities and 
coupled to the TElOl mode for control and readout. 


with standard techniques for four different geometric de¬ 
signs of transmons as shown in Fig. 1. All devices are fab¬ 
ricated on sapphire substrates with identical processes of 
shadow-mask evaporation and lift-ofP^, and therefore are 
assumed to have the same loss tangent for the same type 
of surfaces. All devices have qubit frequency u}/2Tr « 6 
GHz and cavity frequency a;c/27r « 9 GHz. 

Full electromagnetic simulation of surface participa¬ 
tion ratio of transmon qubits faces significant numerical 
challenges due to the large span of length scales. One 
may attempt to model transmon electrodes and any dis¬ 
sipative interface layers as 2D films, and infer pi from a 
surface integral of electric field energy. However, such an 
integral is divergent towards the edge of the films^^. This 
divergence is avoided only when the material thicknesses 
are fully accounted for, as was done in a cross-sectional 



FIG. 2. Illustration of the two-step simulation strat¬ 
egy for computing surface participation ratios, (a) 

Schematic of a transmon qubit with its electrodes color-coded 
into several regions. Grey and yellow represents the perimeter 
and the interior regions of the wide pads of the electrodes re¬ 
spectively, and brown represents the narrow leads. A global 
coarse 3D simulation can accurately determine the electric 
fields across the yellow region, but not in the grey and brown 
region near the edge of the metal, (b) A cross-section view of 
the electrode near a metal edge. The electric field distribution 
within this plane can be computed by a fine 2D simulation, 
(c) A simplified schematic of the region near the Josephson 
junction, which is simulated by a local fine 3D simulation. 
The MA, SA and MS interfaces are defined in (b) and (c). In 
our final account of surface participation ratios, contribution 
from the region within 1 fim from the junction is excluded. 
All drawings are not to scale. 

simulation of transmission line resonators^^’^"^. Without 
a similar translational symmetry, a proper calculation 
of Pi for a transmon qubit generally requires simulation 
of 3D field distribution in mm-sized space with sub-nm 
resolution in critical regions, far exceeding practical com¬ 
putation capacities. 

To overcome the numerical challenges, we employ a 
two-step simulation technique by combining a coarse 3D 
simulation of the entire qubit-cavity system [Fig. 2(a)] 
and fine simulations of representative local regions 
[Fig. 2(b, c)]. A significant part of the surface participa¬ 
tion is associated with regions with highly concentrated 
electric field such as the edges of the electrodes and the 
leads near the junction. We argue that the electric field 
distribution in these regions should have a local scaling 
property independent of the electromagnetic boundary 
conditions far away. These scaling properties can be ob¬ 
tained from simulations of local regions with sub-nm res¬ 
olution and subsequently applied to the global simula¬ 
tion to compute the surface participation ratios^^. We 
assume thicknesses of t = 3 nm and dielectric constants 
of e = 10 for all lossy interfaces for easy comparison with 
a previous simulation of planar resonators^^. Using dif¬ 
ferent assumptions here would rescale the participation 
ratios but not change our conclusions qualitatively. 

Our simulation shows that a significant contribution to 
surface participation arises from the region around the 
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FIG. 3. Surface participation and qnbit lifetime, (a) Spatial distribution of simulated surface participation contribution, 
shown for MS interface of selected transmon designs. Solid curves show cumulative integral of the MS surface participation (pms) 
from the electrode leads as they extend from the junction in distance (horizontal axis), indicating total surface participation 
from the junction to that point. The thick dashed lines indicate the total pms of all features. The dotted lines notionally 
represent contribution from the electrode pads as a whole, (b) Open symbols show measured transmon l/Ti as a function of 
Pms excluding contribution from the sub-micron “near-junction region” (green shaded area in (a)) that most probably contains 
no TLS. Red dashed line is a fit to Eq. (1). The same set of data including the near-junction contribution is plotted as grey 
filled circles, with corresponding fit to Eq. (1) shown as the dotted line, (c) SA, MA surfaces and substrate bulk participation 
ratios {psA, Pma and Pbuik) as a function of pms for transmon devices in this study. Dashed lines are guides to the eye. 


junction leads less than 100 nm away from the junction 
itself [Fig. 3(a)]. This contribution is mostly independent 
of electrode geometry, and can be dominant for devices 
with relatively small surface participation^^. However, 
if surface dielectric dissipation originates from a discrete 
set of TLS with density similar to junction defects’^’^"^^^® 
(~ 1 ^m“^GHz“^), it is most likely that such a small 
volume of macroscopically lossy material contains no res¬ 
onant TLS and thus appear dissipationless. This moti¬ 
vates us to introduce a dimensional cutoff and exclude 
the participation contribution from this near-junction re¬ 
gion. We choose to set this cutoff at a distance of 1 ^m 
from the junction, but any choice on the order of 100 nm 
to 10 /im does not affect the total participation signif¬ 
icantly because the participation contribution from this 
intermediate region of the electrode leads is insignificant 
[Fig. 3(a)]. The resultant total pms from the rest of the 
MS surface is approximately proportional to the mea¬ 
sured 1/Ti for all our devices [Fig. 3(b)]. Similarly, we 
also observe pma and psA proportional to 

The proportionality between qubit decay rate and sur¬ 
face participation ratios strongly suggests surface dielec¬ 
tric loss as the dominant relaxation mechanism for all 
transmons in this study. Based on Eq. (1), any geometry- 
independent dissipation mechanism is expected to induce 
a constant relaxation rate Fg to all our devices. If we 
were to include the near-junction contribution (as noted 
above) in pms, a linear fit of our data to Eq. (1) would 
produce an unphysical negative y-interception [Fig. 3(b)]. 
This reinforces the notion of spatial discreteness of sur¬ 
face loss and the necessity of a cutoff. After implement¬ 
ing the cutoff, we see a very small residual qubit decay 
rate (3 ± 1 ms“^), which can be fully explained by the 
magnitude of quasiparticle dissipation and vortex ac loss 
as we noted previously. Therefore there is no evidence of 
any geometry-independent loss mechanisms, such as from 
the crystalline substrate or the Josephson junction itself, 


that limit transmon lifetimes on the level of Q 10 ^. 
The absence of loss from the junction may be a result of 
the small junction size (0.04 /rm^) so that no resonant 
junction defects are encountered in this study. We also 
note that surface loss mechanisms consistent with our 
observed geometric scaling should not be viewed strictly 
due to impurity or defect-like TLS. Potential alterna¬ 
tive mechanisms closely related to surface electric field 
energy, such as phonon radiation due to surface piezo¬ 
electricity^^’^®, may also be broadly included in the sur¬ 
face dielectric loss in this analysis. 

We cannot determine which of the three surfaces are 
the dominant contributor based on these data alone, 
because all three participation ratios change approxi¬ 
mately in proportion when the qubit geometry is varied 
[Fig. 3(c)]. We can determine a weighted sum of the loss 
tangents of the three surfaces, tan (5 ms + 1.2 tan 1 ) 5 ^ -|- 
O.ltani^MA = (2.6 ± 0.1) x 10“®. To extend our analy¬ 
sis to distinguish different interfaces, one generally needs 
to go beyond a planar layout of transmon electrodes, for 
example by incorporating striplines or microstrips. 

We have further calculated or estimated pms for 
reported planar and 3D transmons from the litera- 
ture®’^®“®^, and Fig. 4 shows the Q factors or Ti’s of 
some of these devices as a function of pms- All data 
points with a single-step aluminum lift-off process sim¬ 
ilar to ours fall near or below the surface-loss line of 
tan(5 = 2.6 x 10 “® (red dashed line), consistent with 
the surface dielectric loss determined in this study. We 
believe similarly-fabricated qubits performing substan¬ 
tially worse than this surface-loss line are limited by other 
mechanisms. Early generation of planar transmons may 
incur losses due to non-equilibrium quasiparticles or lossy 
components of the device package^®, and the 3D “verti¬ 
cal” transmons®^ may be severely limited by conduction 
loss across the cavity seam®®. 

Several recent studies used subtractively-patterned 
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FIG. 4. Transmon lifetime vs. MS surface participa¬ 
tion ratio for selected literature data. Open symbols 
represent various 3D and planar transmons fabricated with 
the single-step aluminum lift-off process similar to this study. 
Filled symbols represent transmons with electrode pads fab¬ 
ricated with etch processes to preserve clean MS surface. The 
vertical axes represent qubit quality factors or the equiva¬ 
lent Ti at Lu/27r = 6 GHz. The bottom axis shows Pms in a 
reversed log scale. The top axis shows the equivalent interdig¬ 
itated capacitor pitch width of a planar qubit for correspond¬ 
ing Pms, a helpful alternative unit of surface participation.^^ 

MBE aluminum® or TiN®® films for transmon electrodes. 
These processes were intended for preserving pristine MS 
interface, and subsequent improvement of Ti suggests MS 
interface may indeed play an important role in the total 
surface loss. We find several data points for these qubits 
(the leftmost filled symbols) above our tan^ = 2.6 x 10“® 
line at relatively high pms, confirming higher surface 
quality than have been measured in this present study. 
However, these surface improvements have not been fully 
translated into the best possible performance for devices 
with lower pms, as indicated by their surface-loss bounds 
(blue and green dashed lines in Fig. 4). It suggests the 
presence of other dissipation channels yet to be fully 
suppressed in these high-material-quality planar qubits. 
These devices also include shadow-mask evaporated junc¬ 
tion leads with lower quality surfaces that can have ap¬ 
preciable surface participation and limit qubit Ti. 

Looking forward, further advance of coherence times 
of superconducting qubits will hinge on a combination of 
improving material surface quality and further reducing 
surface participation ratios. The state-of-the-art planar 
transmons have implemented large-sized planar capac- 
itors®°’®® to reduce surface participation, yielding sub¬ 
stantial gains in qubit lifetimes. One may naively expect 
that millimeter-sized 3D transmons may have smaller pi 
by orders of magnitude and make dielectric loss irrele¬ 
vant. The present study shows this is not the case. Fur¬ 


thermore, our simulations find that merely engineering 
larger and more-separated electrodes will incur signif¬ 
icant Pi from the metal leads required to wire up the 
Josephson junction. Nevertheless, substantial further re¬ 
duction of surface participation in qubits can be achieved 
by more complex three-dimensional designs such as deep- 
etched®® or suspended structures"^®. With no hard limit 
in sight, innovative low-participation designs and im¬ 
proved surface quality, together with modest progress in 
suppressing non-equilibrium quasiparticles, are expected 
to bring another order of magnitude increase in the life¬ 
time of transmon qubits. 
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I. SIMULATION METHODS FOR SURFACE PARTICIPATION 

Participation ratios embody a convenient method to account for dissipative loss in dielectric systems. The partic¬ 
ipation ratio of a certain material or component of the circuit can be calculated by integrating electric field energy 
in an electromagnetic simulation of the exact model of the device. However, from typical adaptive-mesh simulation 
techniques, it is very difficult to produce convergent values of total field energy stored in thin surface layers in a 3D 
qubit-cavity system due to the disparity of length scales. To address this challenge, we introduce a two-step simulation 
technique to calculate participation ratios for three different material surfaces—metal-substrate (MS), substrate-air 
(SA), and metal-air (MA)—for a variety of 3D transmon qubit designs. The results of simulated surface participation 
ratios and measured lifetimes for these qubits (used for Fig. 3 and 4 of the main text) are listed in Table SI. In this 
section, we describe the simulation methods used to obtain these results. 

We use a commercial high-frequency electromagnetic solver (Ansys HFSS) to simulate the entire qubit-cavity system 
on a pm-to-mm scale [Fig. SI (a)], where the aluminum film and surface dielectric layers are modeled as 2D sheets with 
zero thickness. The Josephson junction and the aluminum leads very close to the junction (within 1 pm) are modeled 
as a lumped element. This simulation is carried out at the qubit frequency, and similar to those routinely done for 
black-box quantization of cQED systems^. It provides the overall electric field distribution on a coarse scale (~ pm), 
but does not accurately reflect the highly-concentrated fields at electrode edges or near narrow leads approaching 
the junction that are critical to the total surface participation ratios. To take into account the field distribution in 
these regions and supplement the global simulation, we perform additional local electrostatic simulations (using Ansys 
Maxwell) with sub-nm resolution. 

For convenience, we divide the surface dielectric layers in a transmon qubit into two regions: 1) those associated 
with the large “pads,” metal traces > I pm wide intended to form the external shunting capacitor of the transmon, 
and 2) those associated with the narrow “leads,” metal traces < I pm wide that are used to wire-up the junction 
with the pads. Such definitions are straightforward for MS and MA surfaces in direct contact with the electrodes. 
For the SA surface, we associate SA dielectric within I pm of lead edges with the “leads,” and the remainder with 
the “pads.” 

In this study we mostly vary the geometry of the pads to vary the total surface participation ratio. The observed 
changes in Ti, largely correlated with pad surface participation, highlight their importance. However, in our analysis 
we also explicitly calculate surface participation from the leads. This contribution has not been considered before in 
various implicit applications of surface participation analysis of planar resonators to Josephson-junction qubits^’^. 


TABLE SI. Results of simulated surface participation ratios and measured lifetimes of transmon qubits. Par¬ 
ticipation ratios are multiplied by 10“"^, and calculated by summing over contributions from various regions of the surfaces as 
shown for pms for example. Starred “pi total” excludes contribution from the region within 1 pm from the junction. Measured 
Ti’s are listed for individual devices, with the uncertainty representing one standard deviation of its fluctuation over time. 


Design 

PMS,near ^^MS,mid PMS,far PMS,per PMS,int 

leads leads leads pads pads 

(< 1pm) (1 — 10pm) (> 10pm) perimeter interior 

Pms 

total 

Pms 

total* 

PSA 

total* 

Pma 

total* 

ri(ps) 

A 

2.4 

0.16 

0.01 

0.51 

0.32 

3.4 

0.99 

1.18 

0.11 

75±6, 66±7, 95±8 

B 

2.2 

0.20 

0.19 

0.88 

0.37 

3.8 

1.64 

2.02 

0.19 

34±4, 45±5, 43±3 

Cao 

2.2 

0.17 

0.09 

2.33 

0.82 

5.6 

3.41 

4.01 

0.37 

31±3 

C 20 

2.2 

0.18 

0.09 

2.83 

0.95 

6.3 

4.05 

4.72 

0.44 

25±2.5, 26±2.5 

C 15 

2.2 

0.17 

0.09 

3.37 

1.05 

6.8 

4.59 

5.36 

0.50 

25±2, 16±2, 18±2 

Cio 

2.2 

0.17 

0.09 

4.23 

1.19 

7.9 

5.69 

6.63 

0.64 

19±1.5 

Ce 

2.2 

0.17 

0.09 

6.05 

1.35 

9.9 

7.67 

8.96 

0.89 

11±1 

Ca 

2.2 

0.17 

0.09 

11.1 

1.17 

14.6 

12.4 

14.6 

1.6 

7.5±0.6 

Cl ,5 

2.2 

0.17 

0.09 

19.7 

1.42 

23.6 

21.4 

23.4 

3.2 

5±1 

D 

2.1 

0.19 

1.27 

0.58 

0.36 

4.5 

2.40 

2.82 

0.41 

39±4 
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(a) global 3D simulation 


(b) cross-sectional simulation 



FIG. SI. (a) The global high-frequency 3D simulation includes features in the entire centimeter-scale cavity, such as substrate, 
pads, and most of the leads. Conductors are simulated as perfectly-conducting sheets, and interface layers are omitted. The 
inner (orange) and perimeter (gray) regions, described in Section A, are separated by the dashed green contour. The red window 
shows the placement of the cross-section in (b). (b) This cross-section of the edge of an electrode pad on the substrate (light 
blue) labels the geometry within the 2D electrostatic simulation that supplements the global simulation along pad edges. The 
three interfaces of interest with thickness t are shown in red (MS), purple (MA), and blue (SA), while the two regions of the 
superconductor (thickness h) are shown in orange (@, “interior region”) and gray (“perimeter region”). The perimeter region 
is divided into a cross-hatched region (Q, which fails to converge in the global simulation, and a region @ that is convergent in 
both simulations. The division between @ and (B) occurs at x = xo- When simulated, symmetric boundaries are established to 
represent an interdigitated capacitor (IDC) style device with conductor width w and gap width g {w,g <C Xq). All dimensions 
are not to scale. 


A. Surfaces associated with the electrode pads 

The electrode pads are the large structures of the transmon qubits that determine qubit-cavity coupling and are 
often close to a millimeter in size. Despite their large area, the majority of the electric field energy stored in their 
associated surfaces exists near the edges of the pads. Approximation of the electrode pads and surface dielectric layers 
as 2D sheets, a necessary step in full-scale simulations, results in divergent integrals for total field energy at these 
edges. 

To avoid this divergence, we first divide the electrode pads and the associated MS and MA surfaces into “perimeter 
regions” and “interior regions” [Fig. Sl(a,b)] with their boundary set at a constant distance (xq, typically 1 pm) from 
the edge. (The SA surface can be similarly divided by a contour at a constant distance xq from the outside of the edge. 
The treatment of the SA surface is otherwise analogous to that of MS.) In a global coarse 3D simulation, electric field 
in the interior regions does not have sharp variations, and therefore easily converges to spatial distributions that we 
may immediately record as EMA(a:, y) and EMs(a;, y) at the top and bottom surfaces of the electrode pads respectively. 
We use these field distributions to calculate the surface participation associated with the interior region of the pads 
(denoted by the subscript “int”): 


Pi,int=tJJ^ ^\'Ei{x,y)\‘^dxdy/Utot (SI) 

where i = MS or MA, and C/tot is the total electric field energy in the entire space (dominated by energy in the 
substrate and vacuum). Here we have multiplied the field integral by the assumed thickness of the surface layer, t = 
3 nm, further assuming that the electric field is uniform across that thickness. 

The perimeter regions can be described by a spatial coordinate {x^y,z) as shown in Fig. Sl(a,b), where the y-axis 
winds around the edges of the pads, remaining tangent. We further divide the perimeter regions into two halves. 
Energy in the half adjacent to the edge (0 < a:: < xq/^) fails to converge, regardless of initial mesh parameters, 
following the adaptive mesh refinement process. The other half {xo/2 < x < xq) can be made to converge using 
mesh parameters that are computationally accessible. The key concept to our strategy is to employ a constant ratio. 
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Distance from Edge Along x (^m) 


Perimeter-interior boundary position (pm) 


FIG. S2. (a) Normalized distribution of electric field within a cross section of the MS interface near a metal edge, f{x, —tl2). 
It is calculated from the edge inward along a line that bisects the MS interface [Fig. Sl(b)]. Values are plotted for simulations 
of three sets of boundary conditions, where the first number represents the width of the conducting feature w, and the second 
number is the distance between features, g, as in Fig. Sl(b). (b) A sampling of MS scaling factors Fms [Eq. (S3)] versus 
the defined width of the perimeter region {xo), plotted for various boundary conditions, following the labeling convention of 
(a). The dashed black line marks a;o = 1 pm, the most common choice in our practice. This value strikes a balance between 
boundary condition insensitivity and computational ease. 


or “scaling factor” Fi, to convert the integrated field energy in the convergent half into that of the entire perimeter 
regions, so that 


Pi,per = Fit f dx (f ^\'Ei{x,y)\‘^dy/Utot (S2) 

Jxq/2 Jy ^ 

The spatial distribution of electric field in the perimeter region can be written using separation of variables as 
\E{x,y^z)\ = C{y)f{x,z) in the limit of x, z —>■ 0. This is because the electric field near a metal edge should have a 
local scaling property independent of distant electromagnetic boundary conditions. Here f{x, z) describes the edge 
scaling that can be applied to any cross section, independent of y. The actual form of /(x, z) depends on material 
thicknesses and dielectric constants and is difficult to derive analytically. However, we can compute /(x, z) in a 2D 
cross-sectional electrostatic simulation of an electrode pad, which focuses on the metal edge and takes account of 
the actual thicknesses of each material [Fig. Sl(b)]. The reduced dimensionality allows for accurate computation of 
the field inside the surface layer using sub-nm spatial resolution. In this simulation we choose boundary conditions 
representative of the width of the pad (w) and the spatial separation between the opposing electrodes (g). Although 
such a cross-sectional simulation does not accurately reflect the boundary condition in 3D space, as we already noted, 
f{x,z) is independent of the distant boundary conditions as long as x, z ^ As an illustration, /(x, —1/2) is 

shown in Fig. S2(a) for a few very different values of g and w. 

For our devices with electrode pads typically 10 to 500 pm in their smallest dimension, and separations on about 
the same scale, the above edge scaling function /(x, z) is a very good approximation within the perimeter region for 
properly chosen xq. From /{x, z) we can calculate the scaling factor F) based on the ratios of integrated field energy 
within the cross section: 


fo° dxf°^ p{x,z)dz 

Ixo/2dxI-t P(.x,z)dz 


(S3) 


Fma — 


fo° dx p{x, z)dz + dx pjx, z)dz 
Cel2^^\h^' fHx,z)dz 


(S4) 


Scaling factors Fms for various extents of the perimeter region are shown in Fig. S2(b). We limit our method in 
the regime of xq ^ g, w, where Fms is insensitive to the values of g and w. In practice, we use xq = 1 pm for most of 
the pad structure (which are at least 10 pm in width and separation). Inserting these simulated scaling factors into 
Eq. (S2) allows one to arrive at Pi,per- 
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FIG. S3, (a) A local, near-junction simulation is performed to evaluate lead participation ratios. Energy is integrated within 
cross-sections of the MS interface layer (red). Separately, energy measured along the bisector (yellow) is used for “stitching”: 
comparing local and global simulations in convergent areas to establish a common energy scale, (b) In this top-view of the local 
simulation, the three regions of MS dielectric are color-coded in accordance with the conductor shading in (a). The junction 
(black) is explicitly simulated with thickness t. The near and intermediate regions, including adjacent SA dielectric, comprise 
the “accurate” area of the local simulation (dashed line bounding box); the far region is adversely affected by simulation 
boundaries and is ignored. The stitching region (diagonal lines) is located at the boundary of the intermediate and far regions, 
and at the edge of the local “accurate” area. The local simulation excludes some of the far region (gray), (c) In this zoomed- 
in top-view of the global simulation, the intermediate and far regions are defined, but the near area, extending distance d, 
is approximated as a lumped element. The accurate portion of the leads lies near their center (dashed line bounding box); 
energy near edges or the lumped-element area is divergent. The comparison of stitching extent (shaded area) along the bisector 
(yellow) to that of (b) determines the constant scaling between local and global electric fields. Within the stitching region, 
the cross-sectional energy density can be normalized by the bisector energy density to obtain an energy ratio f{x, z) that is 
independent of y, useful for finding pi,far from Eq. (S6). 


B. Surfaces associated with the junction leads 

A schematic of the Josephson junction and the leads is shown in Fig. S3(a), where x-axis and y-axis are defined 
perpendicular and parallel to the leads, respectively. We divide the surfaces associated with the junction leads into 
three regions based on distance from the junction: the near region (|y| < 1 pm), the intermediate region (1 pm 
< \y\ < 10 pm), and the far region (|y| > 10 pm). The surface participation ratios for these regions are denoted 
by Pi, near, Pi,mid and Pi,far respectively (as shown in Table SI), where i = MS, MA or SA. The near region of the 
leads is not explicitly included in the global simulation. The intermediate and far regions are included in the global 
simulation, but the surface integration of field energy does not converge due to the influence of edges. A scaling factor 
solution akin to that in Section A demands x,z ^ g,w, but the lead is too narrow and too close to the junction to 
satisfy this. We use a supplemental local 3D simulation of the junction leads as shown in Fig. S3(a), which includes 
the thicknesses of all materials, to compute the surface participation of all three regions surrounding the leads. 

This high-resolution local simulation is performed by applying an electrostatic voltage potential between the pair 
of leads across the junction. The boundary of the local simulation is set sufficiently far (typically 25 pm) to ensure 
the calculated field distribution 'Eicicix,y, z) in the the near and intermediate regions is not affected by the type of 
boundary condition used. The overall magnitude of electric field in this local simulation is arbitrarily set by the 
imposed voltage, and must be rescaled by a constant C to be consistent with the field scale of the global simulation 
from which C/tot is obtained. 

This constant C can be determined by comparing Eioc(a;, y,z) with the field distribution in the global simulation 
Egbi(a:, y, z) in a selected overlapping region (“stitching extent”) where both simulations are reliable. In particular, 
we choose the stitching extent as the center line of the leads in the 5 pm < |y| < 10 pm region [Fig. S3(b)]. Such a 
choice avoids the numerical imprecision of the global simulation in areas close to the junction or the edges. It also 
avoids any artificial boundary effects of the local simulation by remaining distant from the boundary. We confirmed 
the two simulations show consistent spatial dependence over this stitching extent, Egbi(0, y, 0) oc Eioc(0, y, 0), and the 
constant C is computed from the ratio of the two. 

Surface participation ratios for the near and intermediate regions of the leads can then be immediately calculated 
by integrating Eioc(a;,y, z) over the volume of interest. For example, 

Pi.near = , y, z)\’^dx dy dz/Utot (S5) 

J J Ji,near ^ 
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The surface participation ratios from the near and intermediate regions are expected to be independent of the design 
of the electrodes, and therefore show very little change among all the devices reported in this study (Table SI). 

On the other hand, Eq. S5 does not apply to lead energies in the far region, which is not fully included in the 
local simulation. To calculate Pi,far we adopt a separation-of-variables approach by noting that |Egib(a;, y, z)| = 
|Egib(0, y, 0)|/(a:, z). Here f{x,z) describes the cross-sectional distribution of electric held in dimensionless units 
(normalized by the held magnitude at the center line of the lead) [Fig. S3(a)]. It can be obtained from the local 
simulation of the junction leads discussed above, which also conhrms that f(x,z) is independent of y for y 3> 1 qm. 
Therefore, 


Pi, far 


rVia-T 
110 M-m 


^|Egib(0,y,0)pdy/C7tot 


\f{x, z)\‘^dxdz 


(S6) 


where the second integral ehectively produces a constant factor that converts the electric held at a single point of the 
center line into energy per unit length along y. This factor is equal to 7.5 x 10“^® m^ for the typical lead width of 1 
p.m. 


II. FABRICATION METHODS 

Fabrication of qubits were performed using the Dolan bridge technique^ on 430 um thick c-plane EFG sapphire 
wafers. After cleaning in acetone and methanol, the wafer was spun with a bilayer of e-beam resist consisting of 550 
nm of MMA EL13 and 70 nm of PMMA A3, then baked at 175 °C. A 13 nm aluminum him was then evaporated 
as an anti-charging layer for electron beam lithography. Patterning of the qubit was done on a 100 kV VISTEC 
EBPG 5000-1- e-beam writer. The anti-charging layer was removed with TMAH, and the wafer was subsequently 
developed for 55 seconds in 1:3 MIBKiIPA followed by a 10 second rinse in IPA. The wafer was then loaded into a 
Plassys e-beam evaporation system (MEB550S or UMS 300). After a 40 W Ar/02 3:1 plasma cleaning for 30 seconds, 
without breaking vacuum, a bi-layer of aluminum (20 nm and 60 nm) was deposited using double-angle evaporation. 
In between the two layers, the junction barrier was grown by thermal oxidation using a Ar/02 85%/15% mixture at 
15 Torr for 12 minutes. Finally, the aluminum was capped with another oxide layer grown with the same mixture at 3 
Torr for 10 minutes. After deposition, liftoff was performed in 60 °G NMP for several hours, then rinsed with acetone 
and methanol. Prior to dicing, a layer of photoresist was spun on the wafer to protect the qubits. After dicing in an 
ADT ProVecturs 7100 dicer, the resist was removed by rinsing in acetone and methanol. 


III. QUBIT LIFETIME VS. SA & MA PARTICIPATION RATIOS 

In the main text, we presented the linear relationship between qubit relaxation rates (1/Ti) and the MS surface 
participation ratio (pms)- Since for all our qubit designs, MA, MS and SA surface participation ratios change 
approximately in proportion [Fig. 3(b) of the main text], qubit l/Ti also shows similar linear relationship with psA 
or pma, as shown in Fig. S4. Assuming SA or MA is the only lossy surface, linear fits to the two data sets indicate 
tan^sA = 2.2 x 10“^ or tani5MA = 2.1 x 10“^ respectively. Since any of the three surfaces can be responsible for the 
qubit relaxation, these values, together with tani5MS = 2.6 x 10“^ obtained from Fig. 3(a) of the main text, should be 



FIG. S4. Qubit relaxation rate (l/7i) as a function of (a) substrate-air (SA), and (b) metal-air (MA) surface participation 
ratios, excluding contribution from the near region of the leads. Dashed lines are linear fits to the data. 
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considered upper bounds for these dielectric surfaces. Furthermore, since a combination of the three has to explain the 
strongly-correlated changes in 1/Ti, we conclude tan^MS + 1-2 tan^sA + 0.1 tan^MA = (2.6 ± 0.1) x 10“^. Linear fits 
in Fig. S4(a,b) also give residual (geometry-independent) relaxation rates of 3± 1 ms“^ and 4± 1 ms“^ respectively, 
consistent with the value obtained from Fig. 3(a). 


IV. ESTIMATE OF SURFACE PARTICIPATION IN PLANAR QUBITS 

In the main text, we placed a number of reported transmon qubits from literature on a diagram of Ti versus pus 
(Fig. 4) for comparison with devices in this study. For most planar transmons included in the figure, we do not have 
complete knowledge of the geometric parameters related to all aspects of their design {e.g. junction leads, coupling to 
the ground plane, device package, etc.) to perform a full-scale simulation. However, since surface participation ratios 
for most planar transmons are dominated by capacitor pads with approximate translational symmetry [i.e. having 
a longitudinal dimension much larger than the lateral dimension), we can estimate their pms from cross-sectional 
simulations alone, similar to the previous work on CPW resonators®’®. Horizontal error bars of ±15% represents 
uncertainties that can be caused by variations in parameters not captured in such a simulation. The planar capacitors 
that we have simulated fall into three styles: interdigitated capacitor (100)®’”^“®, coplanar waveguide (CPW)^’^® and 
coplanar capacitor (CPC)^^. All three styles can be simulated in settings similar to Fig. Sl(b) with different choices 
of boundary conditions. 

In the context of translation-symmetric planar structures, the surface participation ratios are predominantly con¬ 
trolled by the width (w) of the capacitor electrodes and the gap (g) between them. Assuming w and g are varied 
in proportion, the surface participation ratios are approximately inversely proportional to ic or g. (More rigorously. 
Pi oc ln()^)/w, where t* is related to the thicknesses of the metal film (h) and the surface dielectric layers {t). For 
MS interface and for h = 80 nm, t = 3 nm, « 8 nm.) Therefore, it is convenient to express (inverse) surface 
participation ratios in the form of an effective length scale, Weff- We define Weff, or “effective IDC pitch width”, 
of a qubit under study as the width (w) of an IDC structure (with g=w) with identical metal-substrate participa¬ 
tion ratio. This effective width has been used as the top axis in Fig. 4 of the main text, whose relationship to pms 
is calibrated through cross-sectional simulations of IDC structures. We also find the surface participation of CPW 
and CPC structures with g = w are equivalent to IDC structures with WeS ~ l-8g in both cases. An advantage of 
using WeS is that surface participation ratios across different devices can be compared without assuming hypothetical 
thicknesses and dielectric constants of the surface dielectric layers. 

The uncertainties of Ti for these qubits reported from other institutions are based on the stated uncertainties, 
provided sample statistics, or the variations as a function of frequency (for frequency-tunable qubits). 
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